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We study the asymptotic scaling properties of domain wall networks with three different tensions 
in various cosmological epochs. We discuss the conditions under which a scale-invariant evolution of 
the network (which is well established for simpler walls) still applies, and also consider the limiting 
case where defects are locally planar and the curvature is concentrated in the junctions. We present 
detailed quantitative predictions for scaling densities in various contexts, which should be testable 
by means of future high-resolution numerical simulations. 


I. INTRODUCTION 

Phase transitions in the early universe have a num¬ 
ber of inevitable consequences, including the formation 
of topological defects via the Kibble mechanism 0,1- 
The more interesting among these are cosmic strings, for 
which recent developments in fundamental theory have 
led to a resurgence of interest. Specifically, it was noticed 
that the production of cosmic (super) strings is much 
more generic than previously thought. Recent reviews of 
this topic 0,1 provide additional details of these devel¬ 
opments. The potential observability of these strings can 
open a window into physics at the highest energy scales. 

These objects are much more complex than ordinary 
field theory strings. Studying their evolution is a chal¬ 
lenging problem involving physics from different energy 
scales, which is often not fully understood. Quantita¬ 
tive analytic models, even for the simplest strings 01, 
necessarily rely on averaging over some network configu¬ 
rations, which require the introduction of phenomeno¬ 
logical parameters whose values can only be inferred 
by using numerical calibration (M3- Various types of 
well-motivated string networks remain whose cosmolog¬ 
ical evolution is relatively unexplored, including cosmic 
strings carrying charges or currents |ll| (which arise nat¬ 
urally in SUSY models) and models with junctions con¬ 
necting string segments with different tensions, though 
there has been recent progress in the latter [TS], [l3| and 
some relevant numerical studies have been carried out for 
domain walls [TdO ]. 

A more modest but useful approach is to start by 
studying simple toy models, where it may be easier to 
understand the various dynamical mechanisms at play, 
and then using this knowledge for tackling more compli¬ 
cated (but also more realistic) scenarios. Domain walls, 
being the simplest defect (in the sense that they can be 
described by a single scalar field) provide a particularly 
useful testbed, despite being tightly constrained by obser¬ 
vations [3, [l3— though even this case has its challenges 
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and requires large computing resources EMzl. 

Here we study the asymptotic scaling properties of do¬ 
main wall networks with multiple tensions, by generaliz¬ 
ing an analytic model for the simplest walls. Our main 
interest is on the conditions under which a scale-invariant 
evolution of the network (whose presence for the simplest 
walls is well established, both on analytic and numerical 
grounds) still applies, but we will also consider a physi¬ 
cally interesting low-curvature limit. The goal is to use 
this information to shed light on the analogous case for 
cosmic strings, in which case the issue of the generality of 
the scale-invariant asymptotic evolution is still not fully 
clear, although several analytic [Hj, [3,113 and numerical 
efforts [13, 120 I 1221 1 suggest that it applies in a wide range 
of circumstances. 


II. ANALYTIC MODEL 

An analytic model of a defect network starts with the 
microscopic equations of motion (the Nambu-Goto equa¬ 
tions, in the case of strings) and through a suitable aver¬ 
aging process arrives at ’thermodynamic’ evolution equa¬ 
tions. These must include phenomenological terms to ac¬ 
count for defect interactions and energy losses, and for 
their calibration one must resort to numerical simula¬ 
tions. For t/(l) cosmic strings, one obtains the velocity- 
dependent one-scale (VOS) model 00, which has been 
thoroughly tested against simulations. Analogous mod¬ 
els have also been obtained for monopoles [23[ and for 
domain walls [2dj . In what follows, we provide a con¬ 
cise derivation of the analytic model for domain walls, 
summarizing the results of (24| . 

Let us start by considering a flat Friedmann- 
Robertson-Walker (FRW) universe containing a network 
of non-interacting infinite planar domain walls, oriented 
along some fixed direction and all with the same value 
of the velocity v. The momentum per unit comoving 
area of the walls is proportional to a -1 (where a is the 
cosmological scale factor) so that we have V'y oc a~ 3 . 
Differentiating and noting that the Hubble parameter is 
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H = (1 /a)(da/dt), we have 

4 ~ + 317(1 — v 2 )v = 0 . (1) 

at 

In the absence of interactions between walls, the aver¬ 
age number of domain walls in a fixed comoving volume 
should be conserved so that p oc 7 a -1 or equivalently 

^ + H(1 ± 3v 2 )p = 0 , ( 2 ) 

where p is the average energy density in domain walls 
and we have again used Eqn. © to obtain ©• 

Now, instead of assuming the domain walls to be in¬ 
finite and planar we may define a characteristic length 
scale, 


overall energy losses and can account for the limiting case 
where there are no losses (by setting c —> 0 ). 

On the other hand, we also assumed planar infinite 
walls. Realistic wall networks will of course be curved, 
and their curvature will be responsible for an accelera¬ 
tion term which also needs to be taken into account. The 
acceleration should be inversely proportional to the cur¬ 
vature radius, and in the context of a one-scale model the 
curvature radius is identified with L itself. To account for 
the fact that this identification, although approximately 
valid, need not hold exactly, a further dimensionless cur¬ 
vature parameter k (expected to be of order unity) is also 
introduced. The velocity equation is thus corrected to 

!-»■)■ (7) 


L -\' (3) 

which is directly related to the average distance between 
adjacent walls, with a being the domain wall mass per 
unit area. In cosmologically relaistic domain wall net¬ 
works this characteristic length scale L will be approx¬ 
imately equal to the curvature scale of the walls, and 
the walls will have a non-zero probability of crossing and 
interacting with each other. Moreover, the value of the 
velocity will vary along the walls; nevertheless, in the ab¬ 
sence of interactions Eqn. © would remain valid, with 
v being taken as the RMS velocity of the strings. 

We therefore need to add further terms to Eqns. © 
and © in the context of realistic models. Let us con¬ 
sider the latter equation first. The probability of a wall 
element of characteristic size L encountering another seg¬ 
ment of the same size within a time dt is proportional to 
vdt/L [I]. We thus expect an energy loss rate 


loss 


= ~ C L P ’ 


( 4 ) 


where we have introduced a dimensionless proportional¬ 
ity factor, c, which we may expect to be a constant. (This 
expectation is in agreement with numerical simulation re¬ 
sults.) We may therefore add this term to the right hand 
side of ©, which will account for the energy lost from the 
wall network due to the production of wall blobs—also 
called ’vacuum bags’. This mechanism is analogous to 
loop production in cosmic string networks, although the 
vacuum bags will typically dacay much faster than the 
string loops }l5| ] . The evolution equation now becomes 

— + H(l + 3v )p= -c-p = - p, (5) 

dt La 

or, writing it in terms of the length scale L 1 

= (1 + 3 v 2 )HL + cv ( 6 ) 

dt 


This assumes that walls do intercommute when they in¬ 
teract, but the phenomenological parameter c models the 


Note that physically the evolution equation for L (or the 
density) is an energy conservation equation, while that 
for v is analogous to Newton’s second law. 

The above derivation can be analogously done for 
strings, in which case on recovers the evolution equations 
that can be independently obtained startng from the mi¬ 
croscopic equations of motion derived from the Goto- 
Nambu action @,@]- This, together with comarisons with 
numerical simulations, lends further support to the va¬ 
lidity of these evolution equations. We emphasize that 
this is a one-scale model , meaning that one assumes that 
there is a single relevant lengthscale in the problem. 
In addition to the characteristic lengthscale L defined 
above (which is essentially a parametrization of the en¬ 
ergy density in the wall network) one can define a cor¬ 
relation length £ and a curvature radius R, for exam¬ 
ple. In the context of this model, one is assuming that 
L = £ = R. (In the case of cosmic strings, high-resolution 
Goto-Nambu numerical simulations have shown that this 
is a good approximation [7|.) 

The model therefore has two free parameters, c and 
k: the former quantifies (fractional) energy losses by the 
network, while the latter quantifies the curvature-related 
forces acting on the walls. To a first approximation these 
are expected to be constant, and this has been confirmed 
in recent high-resolution numerical simulations Ema, 
which found 

c = 0.34 ±0.16, k = 0.98 ±0.07. ( 8 ) 

Effectively, these provide a numerical calibration for the 
analytic model. 

Neglecting the effect of the wall energy density on the 
background (specifically, on the Friedmann equations), 
which is the usually the case when one carries out high- 
resolution field theory simulations of these networks, one 
can show that the attractor solution to the evolution 
equations m is a linear scaling solution 

L = et , v — const . (9) 

Assuming that the scale factor behaves as a oc t A , the lin- 
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ear scaling parameters have the following detailed form: 


k(k + c) 
3A(1 - A) ’ 


2 1 — A k 

v 2 = -. 

3A k + c 


( 10 ) 


In practice, when carrying out numerical simulations one 
looks for the best fit to the power laws A/V oc p w oc 
l/£ c oc r^ and "/v oc rf, where £ c is the comoving corre¬ 
lation length (the numerical counterpart of the physical 
lengthscale L), A/V is the fraction of the simulation box 
containing walls, and 7 is the Lorentz factor. For a scale- 
invariant behavior, we should have p = — 1 and v = 0. 


III. MULTIPLE TENSIONS 

We now consider a model with domain walls of 3 dif¬ 
ferent types, with tensions (ie, wall masses per unit area) 
oq, (T2, 03. The purpose of this model is to describe the 
interaction between the different types of defects, in order 
to determine if this kind of networks can reach a scaling 
regime, as suggested in [TsL ITol . Additional motivation 
for working with three types of walls is also provided 
by [13, who find that in the case of cosmic strings the 
three lightest components dominate the dynamics, with 
the heavier ones being suppressed. 

In describing a network with 3 types of walls, one might 
start by using a set of evolution equations for L (or p) 
and v for each of the wall types. In fact, for many pur¬ 
poses one can simply use a single, averaged velocity for 
the entire network; we will comment further on this as¬ 
sumption shortly. More importantly, one must have fur¬ 
ther phenomenological terms describing the interaction— 
specifically, the energy transfer—between the three types 
of walls. This will also lead to the appearance of addi¬ 
tional parameters in the model. 

For simplicity, in this model we consider that from the 
3 types of walls, types 1 and 2 behave as standard wall 
networks, with interactions between defects of the same 
type being modelled as in the standard case. But, when 
walls of type 1 collide with those of type 2 (or vice versa), 
they produce a segment of a new type of wall, type 3 
(with tension (73). This behaviour requires the introduc¬ 
tion of an additional term, modelling the interaction be¬ 
tween walls of types 1 and 2. 

We can now model this system by generalizing the one- 
scale model description. Specifically, using the results 
of the previous section (and in particulat the fact that 
Pi = Cj/Lj), the evolution equations for the densities 
of the three types of walls and for the average network 
velocity become 

dpi /, 1 , cv 2 1 • n 1 ■, 

— = -(1 + 3 v)H Pi - p»--p 1,2—>3 (11) 

at di Z 

for z = 1,2; 


and 

£=“-*’)(£- “■)• (13) 

These are therefore the usual evolution equations for 
each type of wall (through written in terms of the density 
rather than correlation length, cf. the previous section), 
except that a common velocity has been assumed for all 
of them (more on this below) and a phenomenological 
interaction term between the different types was added, 
which in our case describes the production of wall seg¬ 
ments of type-3 as a result of the collision of type-1 and 
type-2 walls. All that remains to be done is therefore to 
obtain the specific form of this term. 

Note that above we have assumed that type-1 and 
type-2 networks have equal energy losses as the result 
of the interactions. This assumption is reasonable for 
networks with comparable tensions (which is likely to be 
the realistic case). Its validity is less clear in the case 
where one tension is much larger than the other. In that 
case, if both networks lose comparable lengths, then the 
lighter one will clearly lose less rest energy. On the other 
hand, the lighter network should be moving faster, im¬ 
plying greater losses from kinetic energy. The relative 
contributions of the two effects may well depend on both 
energetic and topological considarations, and therefore 
be mo del-dependent. Such a detailed analysis is beyond 
the scope of our work, but we do caution the reader that 
this equipartition assumption should be tested with nu¬ 
merical simulations. 

We will need to introduce a new dimensionless phe¬ 
nomenological parameter, which we’ll denote d: an effi¬ 
ciency parameter controlling the rate of energy transfer 
into type 3 walls, much in the same way the c parameter 
does for loops (self interaction between walls of the same 
type) 0. Thysic.ally this term should be symmetric in L\ 
and L 2 (or pi and P 2 ) and grow with both of these den- 
sitities, since higher densities of types I and 2 will lead 
to more interactions, as will a higher average velocity of 
the network. For the same reason, it should decrease as 
one increases the density of walls of type 3, P3. For the 
case of strings, this energy loss term was shown in [19| to 
have the form 


p oc E 3 I{L\L 2 ) (14) 


with E 3 being the energy (of type-3 defects) produced 
in a volume L 3 per unit time. On dimensional grounds 
this behaviour should apply for walls, the only salient 
difference being that for the one-scale model for strings 
the correlation length, density and string tension p are 
related by pstring = p/L 2 tring while for walls we have 
p = a/L. Using the latter definition we can therefore 
write the energy transfer term as 


Pl,2->3 — 


dv 


gj P1P2 

of erf p\ 


(15) 


dp 3 

dt 


-(1 + 2>v 2 )Hp 3 - -p\ + p i i2 —> 3 

0.3 


( 12 ) 


The only caveat is that the phenomenological value for 
the parameter d inferred from string simulations need not 
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hold for walls, and in what follows we will therefore treat 
d as a free parameter, though with the expectation that 
it be of order or slightly less than unity. We note, in 
passing, that the values inferred from simulations for the 
model parameters c and k are also different for strings 
and walls. 

At this point it’s useful to summarize the underly¬ 
ing assumptions in our model for multi-tension walls, 
and to compare them with those in existing models 
for multi-tension strings. Indeed, our assumptions do 
differ slightly from previous works for strings Em. 
While E| assumes a single correlation length and ve¬ 
locity for the entire network (allowing for different den¬ 
sities for each string type), El assumes different corre¬ 
lation lengths and velocities for each type (and relates 
correlation lengths and densities through a Brownian as¬ 
sumption). Our assumptions are therefore intermediate 
between these two. 

It must also be emphasized that the terms describ¬ 
ing the interactions between the three wall types satisfy 
energy conservation. This assumption (not explicitly en¬ 
forced by 0 ) is a conservative approach when it comes 
to scaling. For simplicity we only consider energy trans¬ 
fers from two types of walls to the third one, but we have 
checked that our general results would still hold if all 
possible energy transfers were allowed. In other words, 
the broad qualitative features of the solutions (whether 
or not there is scaling, for example) would remain un¬ 
changed, although the exact values of the scaling den¬ 
sities and velocities would of course depend on the nu¬ 
merical values of the parameters. Since the former are 
our main concern here, we have chosen to concentrate on 
the simpler case: discussing a more general one would 
mean added mathematical complexity but little addi¬ 
tional physical insight. 

The use of a single average velocity for the network is 
a good approximation for the velocity of any of its com¬ 
ponents. This is expected from previous work, but we 
have verified that using different velocities for each type 
of wall yields the same qualitative results as the ones we’ll 
describe. A simple comparison of the two approaches is 
shown in Fig. [T] (with all other model parameters being 
the same in both cases). Asymptotically all three compo¬ 
nents reach a common velocity, although the timescale in 
which they do so is different for different wall types. For 
the purposes of the current work, which is concerned with 
asymptotic scaling properties, an average velocity for the 
network is an excellent approximation. (One would only 
need to go beyond this if calculating detailed astrophys- 
ical observables for comparison with experiments.) 

Along these lines, note that assuming that the phe¬ 
nomenological energy loss parameter c is the same for 
all three networks is a reasonable assumption. On the 
other hand, if we were to allow for multiple velocities 
there would be no reason for assuming equal curvature 
parameters k for each string type (so allowing for multi¬ 
ple velocities would greatly increase the number of free 
parameters). 


V 



FIG. 1. Example of the evolution of the individual velocities 
of the walls of type 1 (blue), 2 (red), and 3 (yellow). The 
dashed line represents the average velocity of the network, 
showing it to be a good approximation. 


Having the evolution equations for the network, we can 
proceed to finding its scaling solutions, if they exist. We 
look for solutions of the type 


h = — = d t a 

L 2 = — = £2 1 0 

CO I 

b 1 

II 

CO 

= e 3 t 7 

Pi 

P2 

P3 

(16) 


v = v 0 t s . 


(17) 


The result of this analysis is that two branches of so¬ 
lutions are allowed, with the parameter distinguishing 
them being the curvature parameter k. We now discuss 
them in turn. 

The first type of solution corresponds to k ^ 0, and is 
characterized by 

a =/? = 7 = 1, <5 = 0, (18) 

that is, all walls scale linearly with time, while the veloc¬ 
ity goes to a constant value vo, which can be determined 
based on the other parameters of the system. The wall 
curvature balances the damping caused by the Hubble 
expansion and ensures that the network will reach a con¬ 
stant relativistic velocity. 

The effect of the model parameters on the scaling so¬ 
lution is easy to infer. A larger <73 will make these walls 
more difficult to produce. An increase in the curvature k 
will cause the resulting network to have a higher veloc¬ 
ity vo, and will also make the resulting densities smaller, 
since the increase in velocity will increase the probability 
that loops will form, dissipating energy. A faster rate 
of expansion will damp the system more, decreasing the 
velocity, and having the inverse effect on its densities. 

Regarding energy losses, the value of c determines 
the effectiveness and likelihood of loops being produced, 
therefore its increase will decrease the overall density, and 
will even out the amount of the three types of walls. Fi¬ 
nally, the effectiveness of the interaction between walls 
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of types 1 and 2 to form type 3 walls, described by d, 
affects the density of type 3 walls in relation to the other 
two types. 

The other type of solution, with k = 0, represents in¬ 
stead a system in which the domain walls are conformally 
stretched, 

a = /3 = 7 = A, S = — 3A; (19) 

in this scenario, the velocity decays, so the defects are 
effectively frozen by the expansion. 

At a phenomenological level, this second solution cor¬ 
responds to the case where the defects are locally pla¬ 
nar and curvature is effectively confined to the junctions. 
In this ‘frustrated’ case the walls are (locally) straight, 
therefore there is no mechanism for acceleration through 
curvature, and the walls will be frozen by the cosmologi¬ 
cal expansion. Other than this, the effects of the different 
parameters will be similar to the previous case. 


IV. EXPLORING PARAMETER SPACE 

It is interesting to specifically discuss the behavior of 
the relative densities of the different types of walls. We 
will assume o-i = nay, fix c and k to their standard values, 
and plot the relative densities as a function of the free 
parameters n and d. The plots presented in Figure [2] are 
for A = 0.7, but one can verify that variations in the 
expansion parameter do not change the overall form of 
this dependence. The value of 0-3 was also fixed to lOOcri. 

The results are consistent with intuitive expectations. 
In the first panel the walls with higher tension dominate 
the network, particularly for higher values of d (when 
the rate of production of type-3 walls is higher). In the 
bottom panel we can see the tendency for type-3 walls 
to dominate the network when their tension is higher 
(though note that the network would take a very long 
time to reach this asymptotic state). For n « 50 there 
starts to be a balance between types 2 and 3, while the 
ones with lower tension (type 1) tend to disappear from 
the network. 

To further quantify how the choice of the parameters 
of the system influences the resulting scaling solution, 
we can obtain contour plots of the value of the functions 
Ei = Li/t and v, as function of the parameters k and 
c, for c, k €E] 0 , 1]. For simplicity we assume the limit 
ay = (J2 = 03 here. Figure [3] summarizes this result 
for the linear scaling solution branch; for these plots we 
have assumed a radiation-dominated era (A = 0.5) and 
an interaction parameter d = 0.1. 

This analysis confirms the results of the previous sec¬ 
tion and allows us to check the range of physically accept¬ 
able values for the different parameters. For example, re¬ 
jecting solutions where the network reaches vo 1, the 
region in the bottom-right of the top panel of Figure [3] is 


excluded, narrowing the possible values for the respective 
parameters. 

The k ^ 0 and k = 0 branches can be compared with 
the relevant free parameters being c and d ; this can be 
found in Figure[JJ We can identify the behavior described 
previously, as with the increase in d the quantity L 3 /L 1 
decreases for both cases, and c has an effect that is much 
more noticeable in the case with wall curvature, since the 
loop term is proportional to the velocity. 

Finally we fix c and k to the values inferred from nu¬ 
merical simulations and generically define <72 = nay and 
03 = may, and plot dimensionless density ratios for the 
network as a function of these new parameters, fixing 
some reasonable d. In Figurc[5]we have done that, fixing 
d = 0.1 and A = 0.5 (radiation epoch). The behavior is 
again consistent with intuitive predictions, since in the 
first panel the ratio is almost not affected by the param¬ 
eter m (except for high enough values of d, as expected). 
In the second panel, it is obvious the existence of a line 
n = m, which separates the regions of dominance of walls 
of types 2 and 3. 

V. CONCLUSIONS 

We have explored a simple model for the evolution of 
a multi-tension domain wall network, containing three 
types of interacting walls. By studying the asymptotic 
scaling properties of these networks in generic expand¬ 
ing universes we have confirmed two classes of solutions: 
one corresponds to the well-known linear scaling solution, 
while the other corresponds to the case where the curva¬ 
ture is concentrated in the wall junctions. Our results 
complement those for cosmic strings 0,0, though we 
emphasize that the way we model the network is slightly 
different and aims to model the properties that are more 
easily measurable in field theory numerical simulations. 

Our results support the notion that linear scaling is a 
generic attractor solution for defect networks in a broad 
range of scenarios. On the other hand, we have also con¬ 
firmed that concentrating the defect curvature on the 
defect junctions would be a possible way to obtain a 
conformally stretched network. While such ’frustration’ 
scenarios 0, 0 are quite tightly constrained 0 , this 
highlights a possible physical mechanism that would in 
principle lead to it. Further exploring this scenario will 
require high-resolution numerical simulations, which we 
leave for subsequent work. 
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FIG. 2. Contour plots of the density ratios of walls of types 1 and 2 (top panel) and 3 and 2 (bottom panel), as a function of 
n = <72/cri (horizontal axis) and the interaction term d (vertical axis). 



FIG. 3. Contour plots of the average velocity of the network (top panel) and of er == L\/t (bottom panel) as a function of 
the model parameters k (horizontal axis) and c (vertical axis). 
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FIG. 4. Contour plots of the ratio of the network correlation lengths L3/L1 as a function of c (vertical axis) and d (horizontal 
axis) for case k = 0 (top panel) and k = 1/2 (bottom panel). 



i> » * s a> mi 
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parameters m (vertical axis) and n (horizontal axis). 



(top panel) and ps/p 2 (bottom panel) as a function of the 


[9] C. Ringeval, M. Sakellariadou and F. Bouchet, JCAP 
0702, 023 (2007) astro-ph/0511646 . 

[10] J. J. Blanco-Pillado, K. D. Olum and B. Shlaer, 
Phys. Rev. D 83 (2011) 083514 arXiv:1101.5173 [astro- 
ph.CO]]. 

[11] M. F. Oliveira, A. Avgoustidis, C. J. A. P. Martins, 
Phys.Rev. D85 (2012) 083515. 

[12] J. Urrestilla and A. Vilenkin, JHEP 0802, 037 (2008) 
arXiv:0712.1146 [hep-th]]. 

[13] A. Pourtsidou, A. Avgoustidis, E. .1. Copeland, 
L. Pogosian and D. A. Steer, Phys. Rev. D 83, 063525 
(2011) arXiv:1012.5014 [astro-ph.CO]]. 


[14] P. P. Avelino, C. J. A. P. Martins, J. Menezes, 
R. Menezes, J. C. R. E. Oliveira, Phys.Rev. D78 (2008) 
103508. 

[15] Y. B. Zeldovich, I. Y. Kobzarev, L. B. Okun, Zh. Eksp. 
Teor. Fiz. 67 (1974) 3-11. 

[16] A. M. M. Leite, C. J. A. P. Martins, Phys.Rev. D84 
(2011) 103523. 

[17] A. M. M. Leite, C. J. A. P. Martins, E. P. S. Shellard, 
Phys.Lett. B718 (2013) 740-744. 

[18] S.-H. H. Tye, I. Wasserman, M. Wyman, Phys.Rev. D71 
(2005) 103508. 

[19] A. Avgoustidis, E. Shellard, Phys.Rev. D78 (2008) 

103510. 











[20] P. McGraw, Phys.Rev. D57 (1998) 3317-3339. 

[21] M. Hindmarsh, P. Saffin, JHEP 0608 (2006) 066. 

[22] A. Rajantie, M. Sakellariadou, H. Stoica, JCAP 0711 
(2007) 021. 

[23] C. J. A. P. Martins, A. Achucarro, Phys.Rev. D78 (2008) 
083541. 


[24] P. P. Avelino, C. J. A. P. Martins, J. C. R. E. Oliveira, 
Phys.Rev. D72 (2005) 083506. 

[25] M. Bucher, D. N. Spergel, Phys.Rev. D60 (1999) 043505. 

[26] P. Avelino, C. J. A. P. Martins, J. Menezes, R. Menezes, 
J. C. R. E. Oliveira, Phys.Rev. D73 (2006) 123519. 


